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We report on the observation of enhanced four-wave mixing via crossover res- 
onance in a Doppler broadened cesium vapor. Using a single laser frequency, a 
resonant parametric process in a double-A level configuration is directly excited 
for a specific velocity class. We investigate this process in different saturation 



regimes and demonstrate the possibility of generating intensity correlation and 
ant i- correlation between the probe and conjugate beams. A simple theoretical 
model is developed that accounts qualitatively well to the observed results. 
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1. Introduction 

Four-wave mixing (FWM) is a nonlinear optical process which has been largely investigated 
along the past several decades, both in homogeneously and inhomogeneously broadened 



atomic samples 



QJs 



lH3|]. Since the very early studies of FWM, the possibility to enhance the 
efficiency of this nonlinear process through the use of resonant atomic or molecular media 
has led to the investigation of FWM processes using different atomic level configurations as, 
for instance, two-level 4J, |5j, three-level A [6|, |7|, and four-level double-A schemes |2j. For 
instance, a highly efficient, low intensity, FWM signal has been first observed in a double-A 
level configuration by Hemmer et al. p|. More recently this particular configuration has 
attracted much attention due to the generation of narrow-band photon pairs in atomic en- 
sembles via spontaneous FWM |9J , as well as the efficient generation of pairs of intense light 



beams showing a high degree of intensity squeezing 10|. In those previous experiments 
usually two different laser frequencies are employed to explore the parametric resonance in 
the three- or four-level schemes. One common characteristic of the FWM process using A 
or double-A atomic levels schemes is the possibility of completely eliminating the resonant 
absorption by the medium due to the phenomenon of Electromagnetically Induced Trans- 
parency (EIT) 12|, ll3] , allowing also the control of the medium's refractive index |l4J . Il5| . 



In this work we investigate the FWM process in a sample of thermal cesium atoms 
using a single laser frequency where different atomic level configurations (A and double- 
A) can be directly excited depending where the laser frequency is tuned inside the cesium 
Doppler profile. In particular, we demonstrate that a resonant FWM process in a double-A 
level scheme can be observed via the usual crossover resonance. We also show that the 
FWM spectrum is strongly dependent on the intensity of the pump beams. Moreover, we 
demonstrate that the generated FWM beam can be correlated or anti-correlated to the 
incident probe beam depending mainly on which configuration of atomic levels is excited. It 
is worth mentioning that intensity correlation of two laser beams in a coherently prepared 



atomic medium, under EIT condition, have been observed before 16Hl8|. Although these 
observed intensity correlations and anti-correlations are essentially classical in nature, they 
add new perspectives to the investigation of quantum correlated beams in such atomic 
system, owing to the previous observation of quantum correlated twin photon employing 
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FWM in double-A level schemes |9|. 

The paper is organized as follows: In section II we describe the experimental setup 



and the observed FWM spectra; Section III is devoted to the development of a simple 
theoretical model based on the density matrix formalism to calculate the FWM spectra 
in different saturating regimes. We also present in this section a comparison between the 
calculated and measured FWM spectra. In Section IV we present the experimental results 
for intensity correlations and anti-correlations observed for different laser frequencies, and 
present a qualitative discussion of the observed results. Finally, in section V we present our 
conclusions. 

2. Four wave mixing in cesium vapor: Experiment 

A. Experimental configuration 

The experiment consists in using the so called degenerated backward four wave mixing 
configuration (DBFWM) in which three beams are sent into a vapor and a fourth beam is 
generated by the induced polarization. The geometry of the beams is depicted in Fig. Ufa). 
Two beams, the forward pump (F) and probe (P) beams with opposite circular polarizations, 
are sent copropagating into the cell. Usually the forward pump beam is much stronger than 
the probe beam. A third beam, called backward pump (B) beam, is sent counterpropagating 
with respect to the formers and has a circular polarization opposite to that of the F beam. 
(we chose a+ for F beam and a- for the P and B beams). The beams couple to an effective 
four level atomic system as depicted in Fig. QJb), with the two degenerated ground states 
corresponding to different Zeeman sublevels of the F=3 hyperfine level of Cs( 6S1/2 ) state 
and the excited states corresponding to zeeman sublevels of the hyperfine levels F'=2 and 
F'=3 of Csf 6P3/2) state. All three beams couple to both excited states. The conjugated 
beam (C) is created by the diffraction of beam B into the Zeeman coherence grating created 
by the F and P beams. Due to phase matching and angular momentum conservation, the 
beam C propagates in opposite direction to beam P and with opposite (0+) polarization . 

The experimental setup is depicted in Fig. |2j We use an external cavity diode laser 
emitting around 852 nm, with linewidth below 1 MHz and output power of about 40 mW. 
The main beam is sent through an acousto-optical modulator, with its non-deflected beam 
used as the F beam. The beam deflected in order 1 pass through another acousto-optical 
modulator (with the same RF frequency as the first one) and the new deflected beam in 
order -1 is used as beam P. Such setup allows us to scan the frequency of P (up) around 
the frequency of F (up). The polarization of F is rotated by 90° before it is combined with 
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Fig. 1. a) Geometry of the DBFWM experiment, beams F and P are copropagating and beam 
B is counterpropagating with respect to both. Beam C is generated counterpropagating with 
respect to beam P. b) The effective four level atomic system used in the experiment, with 
the ground levels corresponding to different Zeeman sublevels of the F = 3 hyperfine level of 
Cs( 6S1/2) while the excited states are Zeeman sublevels of the hyperfine F' — 2 and F' — 3 
levels of Csf 6P3/2). F and C have a + circular polarization, while beams B and P have <r_ 
polarization. 

P using a polarizing beam splitter (PBS). The combined beams pass through a A/4 plate 
producing the desired circular polarizations and are then sent into the cell. The beams 
diameter at the cell is around 3 mm. The transmitted beams are converted back to linear 
polarization by a second quarter wave plate and separated by another PBS. After the cube, 
F is retro- reflected to produce the beam B (with ub = wp)- A 50/50 beam splitter is used 
to reflect the beam C that propagates in opposite direction to P. The laser frequency is 
monitored using a saturated absorption setup. 

The cell used in the four wave mixing experiment is a vacuum sealed 12 cm long glass 
cell with diameter of 2.5 cm and filled with Cs at room temperature. We have surrounded 
the cell with a layer of /i-metal in order to reduce spurious effects from residual external 
magnetic fields. 
B. Experimental results 

In Fig. [3] we present typical curves for the detected conjugated-beam intensity (Iq) as a 
function of u>p for fixed values of u)p and low pump power (Pp = 0.1 mW). For all our 
experiments, the pump power will always be given in terms of Pp, with Pb = 0.7 Pp. 
Figure [3](a) was obtained for lof resonant with the F = 3 — > F' = 2 transition (W32). 
Figure[3](c) presents a similar signal obtained for ojp at the crossover resonance between F'=2 
and F'=3 (uj C o)- F° r Fig. [3](b), uip was chosen between these two values, more precisely at 
Uf = u 32 + 2tc x 38 MHz. Note the subnatural linewidth, around 0.8 MHz, of the FWM 
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Fig. 2. Experimental Setup. PBS means polarizing beam splitter and AOM stands for 
acousto-optical modulator. 

signal in the three cases. The height of the peak, however, varies with the position of uf 
throughout the Doppler profile, with two maxima at the positions shown in Figs. [3(a) and 
(c). 

The two peaks in Figs.|3(a) and (c) are associated with different atomic level schemes: the 
peak at transition uf = W32 corresponds mainly to all beams resonant with this transition 
for atoms with zero velocity along the laser beams, and thus to a three-level A scheme; 
the peak at ujf = wco corresponds mainly to the case where P and F are resonant with 
the transition F = 3 — y F' = 2 (F = 3 — y F' = 3) and B is resonant with the transition 
F = 3 — y F' = 3 (F = 3 — > F' = 2) for atoms with a velocity component along the beam's 



axis of v , 



2k 



-2J7), with A F the frequency difference between the two excited 



levels and k the F-beam wavenumber. Thus, for such velocity class, the system corresponds 
to a double- A level scheme. In our experiment, we have then a system that allows us to 
change the effective experimental atomic level configuration from A to a double-A by simply 
tuning finely the laser frequency. 

When we increase the pump power to Pp = 1 mW, the relative intensities of the conju- 
gated signal for different fixed up exhibit drastic changes (Fig. Hj). A strong DBFWM signal 
occurs when uf is between the crossover resonance and U32, around 2n x 38 MHz on the 
blue side of the F = 3 — y F' = 2 transition (see Fig. 11(b)). The signal for this frequency is 
larger than the signal for the ujf = W32 resonance (Fig. 11(a)), and the signal for u>f = ojco 
(Fig. H](c)) is strongly suppressed. The peaks at Figs. H(a) and (c) also present a noticeable 
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Fig. 3. Intensity of the conjugated beam (Ic) when up is scanned around a fixed value of 
ujf for a) ojp = W32, b) u>f = W32 + 27r x 38 MHz, and c) ujp = ojco- The powers used for 
the beams are Pp = lOOfiW, Pp = 0.7 Pp, and Pp = 150fiW. The same vertical scale was 
used in all graphs. 

larger linewidth, indicating stronger power broadening at those conditions. 

In order to investigate in detail this dependence of the DBFWM signal for different laser 
frequencies and power, we measured the peak height (obtained for ujp = up) as a function of 
Up throughout the Doppler profile around the 3 — > 2 and 3 — >• 3 transitions, and for various 
pump powers. The results are shown in Fig. [5J together with the saturated absorption 
signal indicating the value of uf at each position. The condition ujp = uf was monitored by 
measuring a frequency beating in a fast photodiode. From this figure it is possible then to 
capture the broader picture of the behavior highlighted in Figs. |3]and|H We clearly observe 
the two peaks of the DBFWM signal at up = u 32 and up = uco f° r l° w pump power (curve 
(a) for Pp = 100 /iW). In addition, we notice that there are also smaller peaks around the 
two other resonances present in the region scanned by Up: the transition F = 3 — > F' = 3 
(at W33) and the crossover for the transitions F = 3 — > F' = 2 and F = 3 — > F' = 4. We 
understand that the F = 3 — > F' = 3 peak is largely suppressed due to optical pumping 
to F — 4. The smaller height for the F' = 2 — F' = 4 crossover peak should result from a 
combination of optical pumping and the smaller strenght of the transitions to F' = 4. 

When the pump power is increased, the power spectrum changes completely (curves (b)- 
(d) on Figs. [5]). The intensity of the conjugate signal (Ic) at the crossover resonance changes 
from a peak to a dip, with a similar passage from maximum to minimum occuring around 
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Fig. 4. Intensity of conjugated beam (Ic) when the frequency up is scanned around a fixed 
value of uf for a) uf = W32, b) ojf = W32 + 2n x 38 MHz, and c) up = ojco- The powers 
used for the beamas are Pp = 1 mW, Pp = 0.7 Pp, and Pp = 150 fiW. The same vertical 
scale was used in all graphs. 
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Fig. 5. Intensity of the conjugated beam while scanning the laser frequency, keeping up = 
up = u B . The pump power is equal to: a) 100fiW; b) 200[iW; c) 400fiW and d) 2000fiW. 
The probe power for curves (a)-(d) is Pp = 100fiW. The saturated absorption signal used 
to calibrate the frequency scan is shown in (e). 

the frequencies U32 and W33 for higher pump powers. For high pump powers we obseve now 
that the largest signals occur in the frequency region between the crossover and the W32 and 
W33 resonances. 



3. Theoretical model for four wave mixing process 

In order to understand the observed dependence of the DBFWM signal with both frequency 
and power of the pump beams, we developed a theoretical model to calculate the correspond- 
ing DBFWM intensity lineshapes. To simplify the notation, we consider levels |a),|6),|c), 
and \d) as being, respectively, the two degenerated ground states (in our experiment two 
661/2 Zeeman sublevels), and two single Zeeman excited states of levels 6P 3 / 2 F'=2 and 
6P3/2 F'=3. This corresponds to the simplified four-level scheme depicted in Fig. [H already 
used to explain the experimental configuration. According to the beam polarizations, the F 
beam couples level \a) to levels \c) and \d), while the P and B beams couple level \b) to \c) 
and \d). 

1. Formalism 

We write the system Hamiltonian, in the rotating wave approximation, as: 

H = H + V f + Vb + V p , (1) 

with 

h = EH-woi, (2) 

j 
V F = ^ e <<*-^)j:\a)(j\+h.c, (3) 

Vi = fe'« ( -^1 (|&)<c| - |6><d|) + h.c, (4) 

where j = c,d; I = B,P; Tiujj is the energy of level \j) (levels \a) and \b) are considered 
at zero energy); kj is the wavevector of field /; and h.c denotes hermitian conjugated. 
Vlp = — =£- and Hi = ^w 1 {I = B, P) are the modulus of the Rabi frequencies for the various 
transitions. For simplicity, we have considered the electric dipole moment D with the same 
modulus for all transitions, but with an opposite sign for the transitions from \b) to \d) (see 



Ref. 19] for the sign of the different dipole moments). The field frequencies Up and Uj 1 
are written in the moving atomic frame and are related to the frequencies in the laboratory 
frame by the Doppler shift: Up = up — k ■ v = Up and Up = up + k ■ v. 

We have used the density matrix formalism to calculate the polarization terms responsible 
for the generation of the conjugated beam. The evolution of the density matrix is given by 
the Bloch equations: 

^ = -i[tf, P ] + ( ,t.). ( 5) 
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To write the relaxation terms (r.t.) we chose to treat the atomic system as closed with both 
excited states decaying with the same rate T to the ground states. We also simplify the 
problem by considering that the excited state has equal probability to decay to each ground 
state and that the decaying of the coherence p a b is given by a relaxation rate 7, which can 
be related to the atomic transit time in the cross sections of the beams. 

Although in our experiment all the laser beams are cw, the stationary solutions of the 
Bloch equations cannot be obtained simply imposing fy = 0, since there can be some beating 
on density matrix elements involving level \b) because beams B and P can have different 
frequencies (in the atomic frame). To circumvent this problem we have chosen to treat the 
system of equations perturbatively in ftp, such that the beating between ujp and uip can be 
neglected. We though consider the following system of equations: 



A p (o) = _I 
dt h 

dV h 



H + V F + V B ,pW +(r.t.) , (6) 



H Q + V F + V B ,p^> -- V P ,p 



(P) 



1 



■M 



h 

+ (r-t.) , (7) 

where p^ and p( p ) denotes, respectively, density matrix elements in order zero and in order 
one in ftp. 

From Eqs. © and (J7j) we get two coupled systems of equations containing rotating expo- 
nential terms such as e l ^ At ~ A ' r > (where A=F,P,B). To eliminate such terms we introduce 
the transformations [20]: 

n( 0) = y^ -(°)'( a ' fe ) e 4( aw F+ {,w s) i -( a ^+^s)- r 1 (8) 

a,b 

(p) _ sr J p )>( a M 



A? = E 



a,b,c 



ah X 



x i[(a^ t +6^3+cw^)t-(afe P +6fc s +cfc F )-r] /g\ 

where a, b are integer numbers and c = ±1. Introducing such expansion in the optical Bloch 
equations (Eq. ([7])) and looking for stationary solutions, we obtain recurrence relations be- 
tween the different 07- and 07- ' coefficients. The conjugated beam in our four wave 
mixing configuration (see Fig. Hfb)) is generated by the optical coherence terms between 
the excited levels and the ground level \a) with exponential time and space dependence such 



that Ic oc 



<r. 



(^.(1,1.-1) 



+ o, 



(n(l.l,-l) 



. We have solved numerically the derived system of re- 
currence equations and calculated the DBFWM lineshapes by also numerically integrating 
the corresponding optical coherence over the Maxwellian velocity distribution of the atomic 
ensemble. 

2. Comparison with experimental results 

The calculated DBFWM spectrum for the generated conjugated beam is shown in Fig. [6] 

for different values of the Rabi frequency of the pump beams. 



a) 



b) 



_w_ 




Fig. 6. Calculated spectrum for the conjugated beam intensity as a function of the beam F 
frequency (u>f) for various Rabi frequencies of the pump beams Qf = ^b'- a) Qf = 0.01T; 
b) Qf = 0.2r; c) Qf = 0.5T; and d) Qf = l.Or, where T is the excited levels decay rate. 

For low pump powers the calculated conjugated signal, presented on curve (a) of Fig. |6l 
exhibits clear peaks for laser frequencies resonant with transitions to the excited levels (w c ,d) 
and the crossover resonance (uco)- For higher beam intensities the peaks at ujf = oj c ^ are 
deformed and shifted away, while the crossover peak splits in two peaks, with the splitting 
growing with the beam intensities in curves (b) and (c), and being completely resolved in 
(d), in a such way that the intensity spectrum shows a dip around up = ^co- 

By comparing to Fig. [51 we observe then a quite good qualitative agreement of the 
calculated DBFWM lineshapes with the corresponding experimental results. As for the 
differences between theory and experiment, the most important is the asymmetry of the 
experimental curves from one side to the other of the crossover resonance. This indicates that 
the major simplification of our model was not to consider the optical pumping to the F = 4 
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ground state, which should lead to the weaking of the signal involving the F = 3 — > F' = 3 
open transition with respect to the F = 3 — > F' = 2 closed transition. We also did not 
expect to have a quantitative agreement between the experimental and theoretical saturation 
intensities, since we do not consider the whole Zeeman structure in our model. Taking these 
effects into consideration, however, we notice that such good qualitative agreement indicates 
that we have considered in the simplified theory the essential aspects of our experiment, as 
for the effective four-level system and similar magnitudes of the various dipople moments. 
Such theory must contain, then, the key elements to explain the experimentally observed 
behavior. 

The main point to be clarified is the reason for having dips at the different atomic 
resonances, and particularly at ujco-, once the pump power is increased. As for the crossover 
resonance, one possible reason could be some destructive quantum interference, such that 
ajj^ 1 ' - " 1 ' 1 ) ~ ~° r ad ~ > due t° the fact that for uj f = ujco there are two symmetric 
quantum pathways for the coherence to be built (due to the two excited levels). However, 
by examining the individual terms cr^ ^ l '~ and cr ad ' ' issued from our calculations, at 
uf = ojco and for different atomic velocity classes, we ruled out such effect. We understand 
then that the explanation for the dip formation comes from some saturation effect. We 
envision such saturation coming in the following way. At low intensities, the three observed 
peaks in Fig. [6] comes from the matching of the two-photon resonance for the transitions 
induced by F and P with the individual, single-photon resonances of each exciting wave 
with different atomic levels. At high pump powers, however, the single-photon resonances 
at uj c and Ud are dislocated due to large stark shifts induced by the pump beams. Also, 
the different stark shifts induced by F and B destroy the two-photon resonance for the F 
and P fields. In this way, the strongest signal at high pump powers end up occuring for 
nonresonant frequencies, for which the stark shifts are small and at least the two-photon 
resonance for F and P is still present. 

4. Observation of intensity correlation and anti-correlation 

In order to have more insight about the physical mechanisms behind the experimental sit- 
uation of the last section and to gain knowledge on the possibilities of manipulating the 
correlation properties of the fields, we have investigated the occurrence of classical intensity 
correlations between the probe and conjugated beams. We detected then, simultaneously, 
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the transmitted probe intensity and the conjugated beam intensity with a pair of fast pho- 
todetectors connected to an also fast, 400 MHz, oscilloscope (see Fig. [2]). The beams had 
similar powers at the detectors, and both the optical paths and the electronic cables car- 
rying the photocurrents were adjusted to have the same lengths. An example of correlated 
behavior of the conjugated and probe photocurrents is shown in Fig. [3, which displays the 
signals observed on the oscilloscope as a function of time. 
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Fig. 7. Time fluctuation of a) conjugate and b) probe beam intensities. 



In order to characterize such correlations we have calculated the second order correlation 



function between the two signals defined as 

(Slc(t)SI P (t + r)) 



G<-% 



[6I c (t)] 2 )([6I P (t + r)f 



(10) 



where Ip denotes the intensity of the probe beam, and (Q(t)) = // + Q(t)dt/T denotes 
the data average over an integration time T. Such function gives quantitative information 
on the degree of correlation between the two fields, but it has the drawback of mixing 
all frequency responses and thus giving no information about the spectral behavior of the 
system. To acquire such information we have also calculated, thus, its corresponding Fourier 
transform 21]. 

In Fig. Owe show how the crosscorrelation function spectrum is modified for different laser 
frequencies (we keep ujf = wp). The spectra shown are normalized 21], and in Figs. [S^g) and 
(h) we provide, respectively, the mean intensity of the conjugated beam and the saturated 
absorption signal as we scan the laser frequency, indicating the different conditions applied 
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Fig. 8. Normalized correlation spectrum for different laser frequencies (with up = up = ojb)- 
The different frequencies analyzed are show by stars in g) and h), where conjugated intensity 
and saturated absorption signal, respectively, are shown as a function of the scanned laser 
frequency 

in (a)-(f) with respect to these quantities. We see that for most of the laser frequencies 
the crosscorrelation function spectrum has positive values, meaning a correlated behavior. 
For instance, the correlated temporal series shown in Fig. [7] correspond to the spectrum of 
Fig. Wb). A different behavior, however, is observed for Fig. [5(c) for which the two beam 
intensity fluctuations are anti-correlated. Such situation corresponds to the laser frequency 
at the crossover resonance (ujco) between F' = 2 and F' = 3. 

A simple assumption related to propagation effects leading to beams absorption by the 
atomic vapor can qualitatively explain the different correlated behaviors observed for differ- 
ent laser frequencies. First we should note that the parametric FWM process leads always to 
intensity correlations in the fluctuation of the probe and conjugate beams since the photons 
in these modes are generated in pairs. This inherent intensity correlations should compete 
with the intensity correlation or anti-correlations induced by the differential absorption due 
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to different atomic velocity groups. For the laser frequency at the F = 3 — > F' = 2 transi- 
tion (up = CU32) the main contribution to the conjugated beam intensity comes from atoms 
with v z = 0, thus all beams have the same frequency in the atomic frame and couple to 
the same transition. When the laser frequency fluctuates (phase fluctuation), both beams 
P and C (detected in our experiment) simultaneously drive away or get closer to the atomic 
resonance (see Fig. E3(a)). Therefore both beams are equally absorbed by the atomic velocity 
group with v z — which tends to increase their intensity correlation. Now if we consider 
the absorption of these beams by the atomic velocity group v z > 0, around v z — 0, for 
instance, for the situation depicted in ( Fig. H2a)), the conjugate beam will be closer to 
resonance than the probe beam and therefore will be more absorbed by this group of atoms, 
leading to an anti-correlation in the intensity of these beams. However, the absorption of 
these beams by the symmetric velocity group v z < will induce an opposite intensity anti- 
correlation and, as the velocity distribution function is symmetric around v z = 0, the two 
contribution will cancel each other. Therefore, we conclude that for the laser frequency at 
the F = 3 — > F' = 2 transition (up = W32) the intensity fluctuation on the transmission of 
the probe and conjugate beams should be always correlated due to atomic absorption. 

On the other hand, for laser frequencies at the crossover resonance (corresponding to 
the spectrum of Fig. E^c)), the main contribution to the conjugate intensity comes from 
atoms with v z = Vqo — "of 1 - As before, when the laser frequency fluctuates, for instance, for 
the case depicted in [9](b) the absorption of the conjugate and probe beams by this atomic 
velocity group will generate correlation between these two beams. However, differently 
from the previous case, the absorption of the probe and conjugate beams by the atomic 
velocity groups symmetrically placed around Vco, which induces opposite intensity anti- 
correlations, will not cancel each other since the atomic velocity distribution around voo 
is not symmetric, specifically f(v z > 4^) < f(v z < ^f-). Thus, we conclude that one can 
have an overall intensity anti-correlation between the conjugate and the probe beam in this 
situation. 

For laser frequencies between the crossover resonance and a transition to an excited level, 
such as the situation of spectrum [8](b) where up = U32 + 2tt x (38 MHz), the analysis is 
more complicated since the contribution to the generated conjugate intensity comes from 
different velocity groups. Such correlation measurements, thus, provide another evidence for 
the different velocity classes contributing to the signal at various conditions. In this way, it 
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(a) (b) 

F=3 F'=3 



Fig. 9. Field frequencies in the atomic moving frame for different laser frequencies and 
atomic velocity classes, a) ujf ~ W32 and v z = 0; b) u>f ~ ^co an d v z = vco = "^f 1 - 



highlights once again the possibility of choosing from a three- or a four-level configuration 
for the generation of the DBFWM signal at a room temperature vapor. It also demonstrates 
the possibility of manipulating the correlation properties of the generated fields by choosing 
the position of the pump beams in the Doppler profile. 

5. Conclusions 

We have investigate the four-wave mixing process in a cesium thermal vapor and demon- 
strated that this simple Doppler broadened system can provide different atomic level con- 
figurations for observing this nonlinear phenomenon. We also have developed a theoretical 
model that reproduces reasonably well the measured FWM spectra. In particular, the 
possibility of accessing a double -A level scheme with just a single laser frequency can repre- 
sent a considerable simplification for the production of quantum correlated photons, whose 
generation was already demonstrated in similar systems. Towards this direction, we have 
demonstrated experimentally a strong degree of intensity correlation between the incident 
probe beam and the generated conjugated beam. Moreover, we have demonstrated that 
one can switch from correlation to anti-correlation by simply changing the frequency of the 
incident laser. As noted before, correlation properties between two light beams induced in 
a coherent EIT media might find application in magnetometry and atomic clocks and we 
believe our results will certainly broaden this range of applications. 
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